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Q ABSTRACT 

u 

4^ The escape of ionizing radiation from galaxies plays a critical role in the evolution 

of gas in galaxies, and the heating and ionization history of the intergalactic medium. 
2 Here, we present semi-analytic calculations of the escape fraction of ionizing radiation 

for both hydrogen and helium from primordial galaxies, as well as analytic derivations 
I ^1 of these quantities. We consider variations in the galaxy density profile, source type, 

^ location, and spectrum, and gas clumping/distribution factors. For sufficiently hard 

^ first-light sources, the helium ionization fronts closely track or even advance beyond 

that of hydrogen. Key new results in this work include calculations of the escape 
fractions for He I and He II ionizing radiation, and the impact of partial ionization from 
X-rays from early AGN or stellar clusters on the escape fractions from primordial halos. 
When factoring in frequency-dependent effects, we find that X-rays play an important 
^ role in boosting the escape fractions for both hydrogen and helium, but especially for 

^ He II. We briefly discuss the implications of these results for recent observations of 

^ the He ll reionization epoch at low redshifts, as well as the UV data and emission- line 

^ signatures from early galaxies anticipated from future satellite missions. 

1. INTRODUCTION 

The escape fraction of ionizing radiation, /esc, is a key quantity in the feedback from galaxies 
on the intergalactic medium (IGM) at all redshifts, z, and is an important input parameter in 
reionization models and cosmological simulations. It is typically defined as the fraction of hydrogen 
(or when relevant, helium) ionizing photons that escape from a source galaxy to the IGM after 
accounting for recombinations within the galaxy and its halo. Such radiation profoundly affects 
the thermal and ionization history of the IGM and is an important factor in the modeling of many 
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astrophysical problems, ranging from metal ionization states in absorber systems to the reionization 
of the IGM (Loeb &: Barkana 2001). The converse problem - that of trapped radiation in galaxies 
- is also important for understanding the viability of detecting primordial galaxies or Lya emitters 
through emission-line and other signatures (Tumlinson et al. 2001 Rhoads et al.|[2004 ). 



Theoretical studies (] 


3ove et al. 2000; Wood & Loeb 2000; Ricotti & ShuhpOOO Ricotti et al. 


2002 


Gnedin et al. 2008 


Wise &: Cen 2009 Razoumov &: Sommer-Larsen 2010; jYajima et al. 


2011 


Fernandez & Shull||2011) have indicated a wide range of values for /esc for hydrogen, ranging 



from 1-100%, beginning with the first detailed studies of the Galaxy's Reynolds layer by Dove & 
Shull (1994). Owing to the complexity and variety of factors that influence photon escape, any 



single theoretical trend must necessarily be seen as only one aspect of reality. These factors include 
galaxy mass and morphology, galaxy redshift, gas density profiles, composition, source luminosi- 
ties/spectra, source distribution and time evolution within the galaxy, and clump mass/distribution 
factors in the interstellar medium (ISM). Several numerical simulation studies have recently found 
that /esc increases with decreasing galaxy mass, and with galaxy redshift (Gnedin et al.||2008 Shull 



et al.|2012a| . Higher mass galaxies will have more gas to ionize (leading to lower fcsc), but owing to 
their deeper potential wells, they are also more likely to have higher star formation rates and more 
collisionally ionized gas in their ISM (leading overall to higher /esc, even if the sources are embed- 
ded in neutral or molecular gas). Similar complexities arise in redshift studies of fesc- in principle, 
the abundance of lower-mass galaxies and higher star formation rates at high redshifts should lead 
to higher values of the escape fraction. Some authors find higher values of /esc for dwarf galaxies 
and Lya emitters at high redshifts relative to massive galaxies at similar or lower redshifts, but 



these are complicated by model assumptions underlying these rather different systems (Ricotti & 
Shull||2000[ [Wood Loebl|2000t [Gnedin et al.||2008[ ). Direct comparisons with data of low-redshift 
systems will provide key insights into the many convolved factors, but this has been hindered by a 
conspicuous lack of numerical simulations of /esc for low-z systems. In addition, values of /esc that 
exceed 20% in low-mass galaxies are possibly needed to explain the observational constraints on 
completing hydrogen reionization by z ~ 7 ( Shull et al.||2012a ). The results of the many theoretical 
studies, as well as various factors and trends, are summarized comprehensively in Fernandez & 



Shuh (2011) 



Observationally, the constraints for fesc have so far been for H alone, focussing (through various 
techniques) on measuring the escape of Lyman-continuum (LyC) radiation (hi/ > 1 Ryd). Here, 
there has been a gradual approach to some consensus values, and an indication from the past 
few years of data that /esc,H increases with increasing redshift, confirming one of the theoretical 
predictions stated above. Although observations about a decade ago ( Steidel et al.|[200T ) of Lyman 
continuum emission from Lyman-break galaxies (LBGs) at z ~ 3 indicated high values of /esc,H 



> 0.5, more recent analyses of these objects at similar or lower redshifts (Shapley et al. 2006 



Siana et al. 2010 Nestor et al. 2011) have derived lower values of ~ 0.1-0.25, with only ~10% 
of LBGs and Lyman-a emitters exhibiting large fesc values. Note that these /esc,H values are a 
measurement of the relative escape fraction, the ratio of fractions of escaping ionizing to non-ionizing 
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UV photons, as defined by these authors. There are also indications from measurements of the UV 
luminosity density of galaxy samples at 2; ~ 6-8 that /esc,H may exceed 50% (Finkelstein et al. 



2012), in agreement with Shull et al. (2012a) and with simulations finding increasing /esc values 
towards the faint end of the galaxy luminosity function ( Ciardi et al.|20l"2 ). At lower redshifts {z < 



3), observations have indicated /esc.H ~ 0.01-0.1 from high-mass galaxies (Leitherer et al. 1995 



Deharveng et al. 2001; Heckman et al. 2011 Grimes et al. 2009). Although these fractions appear 



significantly lower than the 2; ~ 3 values at first, at least some of these systems (Heckman et al. 



2011) exhibit values of about 30-40% when recast in terms of the relative escape fraction. Again, 



it would be valuable to have specific theoretical constraints for comparison, but the field has had 
very few numerical simulations of /esc at low redshifts. 

We note that this relative escape fraction is arguably more comparable to the results of this 
work (and many others cited above) than the absolute escape fraction. The absolute escape frac- 
tion includes, by definition, the effects of any dust obscuration. However, compared to hydrogen 
photoelectric absorption, dust is usually unimportant for the EUV (LyC). Our models do not at 
present consider the effects of dust, as they deal with primordial gas. By considering the relative 
escape fraction (usually the ratio of escape fractions at 900 A and 1500 A) the effects of dust are 
mitigated, since the obscuration due to dust should be the same at the two wavelengths modulo 
the difference in wavelength-dependent opacitjj^ 

In contrast, there are no data and detailed theoretical studies of the relevant escape fraction 
for radiation capable of ionizing He 11 {hi^ > 4 Ryd), /esc.He- In this paper, we derive the values and 
geometry of the escape of helium ionizing radiation from the first galaxies. This quantity can have 
a profound influence on the thermal and ionization conditions in the intergalactic medium (IGM) 
and be an important input for future IGM and cosmic microwave background (CMB) studies, as 
well as searches for the first stellar clusters. We find that the hardness of the source spectrum 
is an important factor in determining the relative ionization- front (hereafter, "I-front") evolution, 
and we consider both AGN/QSOs as well as metal- free stellar populations which, unlike their low- 
metallicity counterparts, produce helium ionizing radiation. We also consider the role played by 
X-rays in the escape of ionizing radiation from early galaxies, particularly for He. We present 
our model and assumptions in §2, and the role of frequency-dependent factors in the relative 
propagation of H and He I-fronts in §3. We present our results and detailed calculations in §4, 
allowing for variations in galaxy density profile, source spectrum, source location, and clumping 
and gas cloud factors, and we conclude in §5. In the Appendix, we present an analytic derivation 
of /esc,H and /esc,He for the most relevant subset of cases considered in this paper. 



^The difference in opacity between 1500 A and 1000 A may typically be a factor of 1.85 for the standard ISM 
extinction curve. Future theoretical studies of escape fractions should attempt to incorporate the effects of dust and 
predict relative escape fractions directly. See Benson et al. ( 2001 1 for an example of escape fraction calculations which 
do account for the effects of dust. 
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2. BACKGROUND AND MODEL 

Our goal is to compute the emergent spectrunj^ 5oo(i^), of ionizing photons emitted by some 
source (a star cluster or active galactic nucleus; AGN) with intrinsic spectrum Sq{i^) which is 
embedded within the interstellar medium (ISM) of a galaxy and/or a more extended gaseous halo. 
Given the emergent spectrum and the known input spectrum the escape fraction is defined as: 

f — J 1^0 ' /1^ 

J esc — j-oo o t \A ' V ' 



V 



where vq is the frequency corresponding to the ionization potential for whichever species is of 
interest, with hvQ = 13.6 eV, 24.6 eV and 54.4 eV for H I, He i and He ii respectively. We will 
consider the gaseous environment of the source to be a mixture of hydrogen and helium, and 
to have some density distribution nj(x) for species i = H,He, with the density distributions for 
hydrogen and helium differing only by a constant proportionality factor appropriate for primordial 
abundances. 

We consider a source located at position xq within this density distribution. Along a direction n 
from this source we consider the propagation of the input spectrum through the absorbing gas. We 
assume that at each distance, d, from the source the hydrogen and helium reach photoionization 
balance at a temperature Tq. Given the electron density, ne, the equations of photoionization 
balance imply that 

= SHiraHll - XHii) - a^n(ro)ne?7.HaJHii, (2) 

= SHemHe(l " a^Hell " a^Helll) " "He I (^o)?^enHea^He II 

-EHeIinHe2;HeII + "He iH (-^O^'^enHea^He III , (3) 
= SHcimHea^Hell " aHeiIl(^o)"'enHe2;HeIII, (4) 

(5) 

where the densities are evaluated at position xq + dn, and is 

/•CO 

= / S{u)ai{u)dv. (6) 

J UQ 

(2) 

Here, af'iT) is the case-B recombination rate coefficient for ionization state i, ai{u) is the pho- 
toionization cross-section for ionization state i =H I, He i. He II, and Xi is the fraction of each 



^Typically, the emergent spectrum is defined as that at infinite distance from the source, but without including 
radiative transfer through the IGM. In cases where we consider density profiles that extend to infinity (and whose 
density falls only slowly with distance) we will define the emergent spectrum as that at the virial radius of the dark 
matter halo associated with the source. Our justification is that the virial radius marks the approximate boundary 
between galactic and intergalactic environment, and the escape fraction is usually thought of as the fraction of ionizing 
photons emitted by the source which arrive at the IGM. 
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(2) 

atomic species in ionization state i. Defining TZi = T,i/al_^^nQ, with the convention that 



i + 1 



H II for i = H I 
Hell fori = HeI 
He III for i = Hell, 



these equations have solution 



a^Hii 
a^Heii 
a^Heiii 



'^HII 
1 + ^HIl' 



1 + T^Hell + '^Hell'^Helll 

'^HcIl'^Helll 

1 + "^Hell + '^Hell'^Helll 



(7) 

(8) 
(9) 
(10) 



We will assume a constant gas temperature of Tq = 10^ K throughout, and do not explicitly 
solve for the temperature. We do not expect the temperature within the I- fronts to be much below 
10^ K; in fact, photoheated primordial gas has temperatures > 20,000 K, which has even been 
measured in the Lya forest (Becker et al. 2011). For the lO^M© halo at z = 10 that we will 



consider in ^4.4 the virial temperature is ~ 1.1 x 10^ K ( [Donahue k Shull|l99lt [Becker et al.|201l[ ). 



The ionization balance depends on the electron density, which can be found directly from the 
states: 

He = nH3;HII + nHe(a:^Heii + 2xHoIIl)- (11) 

Therefore, we make an initial guess at the ionization states to allow us to compute Ue and then 
iteratively update the ionization fractions using the above equations until we obtain a converged 
solution for rig. 

After finding the ionization states we compute the optical depth to photoionization across a 
small increment in distance 6d using 



poo 

6t = Sd''^^ / niiXi(Ti{v)<lv, 



(12) 



where i' =H for i =H I and i' =He for i =He i or He ii. The step in distance, 6d is chosen to 
ensure 6t <^ 1 and 6d <^ rg where rg is the characteristic length scale of the density profile being 
considered. The spectrum is then updated using 



S{u) S{u)exp{-5T). 



(13) 



We repeat this process until a sufficiently large distance is reached that the spectrum is no longer 
changing significantly or, when considering source embedded in dark matter halos, the IGM (i.e. 
the halo virial radius) is reached. We then compute the escape fractions, /j_esc(n), as defined above. 
The mean escape fraction, averaged over all directions, is then 

1 

4tt 



{fi,e 



/i,csc(n)dil. 



(14) 
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2.1. Inclusion of clouds 



To include clouds (i.e. regions of enhanced density) in our calculations we follow Fernandez & 



Shull (2011). Clouds are defined by three quantities: their radius (rc), their volume filling factor 
(/c), and the "clumping factor" {C, the ratio of density in clouds to non-cloud regions). Note that 
this definition of the clumping factor is different from the usual one of C = (n^)/(n)^. 

Along each direction, n, we construct a cylindrical volume with radius equal to the cloud 
radius such that any cloud whose center lies within this cylinder will intersect the line of sight 
along which the escape fraction is being computed. We find the mean number of clouds in this 
cylinder based on their radius and volume filling factor, {N^) = /c'/rr^(iniax/(47r/3)r^, where dmax 
is the maximum distance from the source to be considered. The actual number of clouds in the 
cylinder is then drawn from a Poisson distribution with this mean. Each cloud is located randomly 
within the cylinder and the section of the line of sight that the cloud intersects is determined. 

In propagating the spectrum through this cloudy medium we now choose the steps in distance 
5d such that cloud boundaries always precisely coincide with a step. Therefore any given step is 
either in a cloud or not in a cloud (i.e. a step is never partly in a cloud and partly not). In cloud 
regions the gas density is increased appropriateljj^ (see Fernandez & Shull 2011). The emergent 



spectrum of photons then accounts for the effects of clouds along the line of sight. 



3. THE EVOLUTION OF HELIUM IONIZATION FRONTS 



We begin by calculating the critical spectral index at which the rate of H and He photoion- 
izations become equal. For a source spectrum whose specific luminosity (erg s~^ Hz~^) goes as 
, we find that the ratio of photoionization rates per atom for He and H, Xc, is: 



Q(HeII)/nHe 
Q(HI)/nH 



4Ryd/h 



duLy/ {hv) 



/iRyd/h dvL^/{hu) 



"-He 



(15) 



The primordial helium mass fraction Ine ~ 0.25 (Peimbert et al. 2007[ Aver et al. 2010), and the 
currently permitted range of Ob/i^ = 0.02255 ± 0.00054 (Komatsu et al. 2011). This leads to a 
range of values for the relative helium fraction by number, y = nHe/?T-H = 0.0789-0.0823. We find 
that Xc is unity when: 



0.08, or ttcrit 



ln(y) 
'ln(4) 



1.82 ±0.03, 



(16) 



^It is possible for clouds to overlap. In such cases we still increase the density only by a factor of C relative to the 
background, not by a factor of 2C for example. This overlap means that the actual volume filling factor, /c, will be 
slightly smaller than the input filling factor, /ci /c = 1 — exp(— /c) on average. For the cases we consider, /c = 0.2, 
so = 0.18, a small difference. 
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Fig. 1. — The specific intensity as a function of energy of a fiducial 10^ Mq metal- free stellar cluster, 
based on [Tumlinson et al.] ( |2004[ ) . Red and blue lines represent 1-100 M© and 10-140 Mq IMFs, 
with solid and dashed lines in each case denoting the values on the zero-age main sequence (ZAMS) 
and at 2 Myr. Note the relative flatness (hardness) of the spectrum between the H I and He ii 
ionizing thresholds at 13.6 eV and 54.4 eV respectively, and the rapid decline of the He ionizing 
flux with time. For comparison, a typical QSO spectrum (arbitrarily normalized) with slope a = 
1.8 and 0.8 in the respective energy ranges 13.6-300 eV and > 300 eV is also shown. 



if we assume y = 0.080 it 0.003. Therefore, in the case of a source with a < 1.82, the helium I-fronts 
may begin to overtake or coincide with the H I-fronts; such a hard spectral index may already be 
seen in z ~ 2.5-3 QSO spectra (Shull et al. 2012b) as we discuss below. For a source spectrum 
with multiple indices, e.g., qi between vi and V2-, and a2 beyond 1^2, equation (15) is unity when: 



,02 

_2_ 
,ai 



X ^ X 4-°2 

02 



0.08. 



(17) 



To extend this estimate to the case of metal-free stars, we now consider three cases involving 



metal- free IMFs in the mass ranges 1-100 Mq (Tumhnson et al. 2003), 10-140 (Tumhnson 



et al.| ( |2004D ; hereafter TVS04), and > 100 Mq ( |Bromm et al.|200lD , all with a Salpeter slope unless 
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noted otherwise. The first represents a roughly present-day initial mass function (IMF), the second 
the best-case IMF that matches the current data on nucleosynthesis and reionization, and the third, 
an IMF consisting of very massive stars (VMSs). Such stars can lead to pair-instability SNe in the 
mass range 140-260 Mq and may dominate in primordial gas that lacks metal coolants. Thus 140 
Mq divides stellar masses into two regimes that have similar radiative but different nucleosynthetic 
yield properties. For the first two of these IMFs, we display the specific intensity of a fiducial 
10^ M0 zero-metal (Pop III) cluster as a function of energy in Figure [l] In each of these cases, two 
curves corresponding to the zero-age main sequence (ZAMS) and at times of 2 Myr are shown; for 
all the curves, we observe the relative hardness of the cluster's radiation between 1 Ryd and 4 Ryd. 



A numerical evaluation of Equation (15) yields a value Xc — 0.24 (0.07) for the ZAMS (2 Myr) 



1-100 Mq IMF, and ~ 0.43 (0.11) for the ZAMS (2 Myr) 10-140 Mq IMF. 

The fact that Xc never exceeds unity implies that the rate of He photoionizations, and therefore 
the He I-front evolution, cannot exceed that of H for these IMFs. This is mostly due to the sharp 
fall off of the intensity beyond the He ionization threshold. For comparison. Figure 1 also shows 
the spectrum of a quasar (arbitrarily normalized) with slope a = 1.8 and 0.8 in the energy ranges 
13.6-300 eV and > 300 eV respectively ( |Venkatesan et al.||200lD . Such a source would clearly meet 



the criteria in Equation (16), and the He I-fronts could "keep up" with those associated with H. 



AGN with sufficiently hard spectral indices are already observed in COS spectra of 2 ~ 0.03-1.4 



QSOs (Shull et al. 2012b) where composite spectral slopes a = 1.4-1.5 are observed at energies 
1.0-1.5 Ryd. This is a harder spectral slope than the Ocrit derived earlier, if it continues to the 4 
Ryd edge for He 11 reionization. 

In the Appendix, we present analytic calculations of the escape fractions for hydrogen and 
helium as a function of source spectrum frequency for a Navarro-Frenk- White (NEW) background 
galaxy profile, and briefly compare those values with our numerical results below. 



4. RESULTS 
4.1. Gaussian Distribution Model 

We begin by considering the simple Gaussian density distribution model of |Dove &: Shull 



(1994). We use this paper (rather than the many others cited in earlier sections) for general 
comparison, as its analytic results are straightforward to reproduce, at least for a pure hydrogen 
case. It also permits easy comparison with our results which include a number of variations in 
galaxy density profiles and source spectra, as well as the specific new scenarios including helium 



and X-rays. Adopting the infinite slab model given in Equation (7) of Dove & Shull (1994), we 
assume that: 

nH(x) = noexp(-[Z//i]V2), (18) 

where uq = 0.312 cm^'^, the Gaussian scale height h = 0.215 kpc, and Z is the vertical distance 
above the midplane of the slab. The density profile for helium is the same, scaled by the cosmic 
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helium to hydrogen ratio. We adopt a blackbody input spectrum with a temperature of 10^'^ K 
normahzed such that the production rate of H I ionizing photons is similar to that of either a 10^ 
or W^Mq cluster of Pop III stars. In this model, the ionizing photon rates (Q) for a 10^ Mq 
Pop III cluster for H I, He i and He ii are respectively, 1.4 x 10^^ s~^, 6.1 x 10^° s~^, and 2.9 x 10^^ 
s~^, with the Q values for the 10^ Mq Pop III cluster simply scaled up by a factor of 10. In this 
model, there are relatively few X-rays present, although we consider scenarios below with AGN 
that do explicitly include them. We have also considered the contribution from high-mass X-ray 
binaries (HMXBs) which may form nearly simultaneously with the Pop III stars. We adopt the 



model of Mirabel et al. (2011) to specify the normalization of the HMXB spectrum. Specifically, 
their Equation (5) gives the luminosity in the 2-10 keV range per unit star formation rate. We 
find that the inclusion of HMXBs makes negligible difference to the escape fractions and so do not 
consider them further in this work. 



For these Q- values, the analytic formula in Equation (14) from Dove & Shull (1994) for the 
critical opening angle, 9c, of the cone of escaping H- ionizing radiation yields values of approximately 
70° and 81° for the 10"^ Mq and 10^ Mq Pop III clusters respectively. This is in excellent agreement 
with the results of the semi-analytic calculation described in ^ as shown in Fig. [2] Specifically, 
there is close agreement with the angles at which /esc begins to drop significantly below unity, i.e., 
the ionizing radiation is being containecQ The black lines in Fig. [2] show the escape fraction for 
H II when we remove helium from our calculation - clearly the presence of helium causes only a 
small change in the escape of hydrogen-ionizing photons in this particular case. 

Repeating this calculation for He 11 and He iii (with their respective recombination coefficients 
and Q-values as stated earher), we find that for He 11 (He iii), 6c ~ 77° (11°) and 84° (63°) for the 
10^ Mq and 10^ Mq Pop III clusters. This assumes that all of the helium (for yue = 0.08) is either 
in He 11 or He iii for these 9c values. Comparing with the curves in Fig. [2| we see that, similar 
to hydrogen, the 9c for each case corresponds to where the curves depart from their maximum 
value, although this is less well-defined for He iii. This perhaps reveals the limitations of extending 



the Dove Sz Shull ( 1994 ) analysis to helium, especially when we do not account for (in this simple 
estimate) the absorption of helium-ionizing photons by hydrogen. We do see that the critical angle 
for escaping He 11 ionizing radiation is less than that for H I, as predicted by the analytic results 
earlier. 

In Figure [3j we show the escape and ionization fractions as a function of radial distance from 
the source, evaluated at the critical angle of 6* ~ 85°. Note that in Figure [3| the He 11 fraction rises 
again at large radii (> 10 kpc). The density drops so rapidly at these large radii that even the 



' Dove & Shull ( 1994 \ use a Stromgren argument to determine whether any photons along a given line of sight 
can escape. As such, their escape fraction as a function of angle is a monotonically declining function, going from 1 
to at the critical angle. In our calculation we account for the non-zero mean free path of photons through the gas. 
This results in a smooth dependence of escape fraction on angle, making a comparison with Dove & Shull ( 1994 1 
non-trivial. Given this caveat, we find the agreement between our results and the critical angle calculation of Dove| 



& Shull ( 1994 1 to be acceptable. 
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weak radiation field that has escaped to this distance is able to ionize the gas. While interesting in 
its own regard, this is mostly irrelevant for the escape fraction calculation since the density (and 
optical depth) are extremely low at these large radii, especially given our constant temperature 
assumption. (At these low densities, cooling rates will be very low, driving the gas to higher 
temperatures at which collisional ionization becomes relevant.) 

Last, we comment on the geometry of ionized helium versus hydrogen bubbles in the context 
of Dove Sz Shull (1994). Beginning with their Equation (12), we can recast it as, /esc,i(^) = 1 — Wi, 



where Wi oc (a^ nij/Qi^ Here, i corresponds to either H or He, and the constant of proportionality 
is related to the product of the gas distribution and the electron density (both of which are the 
same for H and He). The product (a) nj) is already set for each species, and is typically less for 
He II and He III relative to H ii, whereas for typical Pop HI clusters, Qh significantly exceeds Qne- 
Therefore, in principle, Wue > Wn, and /esc.He < /esc.H- Additionally, extending Equation (14) in 



Dove & ShuU (1994), 



cos{9, 



c,HJ 



cos(6', 



cHei 



Wue 



1/3 



(19) 



which will typically have values less than 1, therefore implying that 0c, He < ^c.H- Thus, the escape 
of He-ionizing radiation will occur through narrower opening angles than for H-ionizing radiation. 
These broad conclusions will hold in the absence of any coupling of their ionization equilibria 
through, e.g., X-rays or otherwise. 



4.2. Blackbody source and clumping/cloud factors 

We next consider the effects of adding in clouds to our model. The left panel of Figure |4] shows 



the escape fraction as a function of angle for the 10^ Mq Pop HI cluster model of ^4.1 but now with 
a single realizatiorj^ of a population of 30 pc radius spherical clouds with a volume filling factor of 
20% and a clumping factor of 10. The denser cloud regions stay more neutral and provide greater 
optical depth to the ionizing photons. Therefore, the escape fractions are lower than the case with 
no clouds. 

The right panel of Figure |4] shows the radial variation of escape fractions and ionization frac- 
tions for the same model. Large fluctuations in ionization state are apparent — these are the differ- 
ence between in-cloud and out-of-cloud points. The escape fraction lines show corresponding steps 
as the photons get strongly depleted when they encounter a cloud. 

The left-panel of Figure [5] shows the escape fraction as a function of angle for ten realizations 

^Note that we use the same radial distribution of clouds for each angle, 9, so there are no fluctuations in escape 
fraction with angle. This is unrealistic of course, but this case serves only to illustrate the effects of clouds. Our 
primary interest is in the mean escape fraction averaged over many realizations of the cloud population, for which 
the lack of angular fluctuations in our calculations is irrelevant. 
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Fig. 2. — The escape fraction as a function of angle, 6, in the Dove & Shull (1994) Gaussian model. 
Thick and thin lines correspond to W^Mq and IO^Mq Pop III clusters respectively. Yellow (solid), 
blue (dashed) and green (dot-dashed) correspond to H I, He i and He il respectively. The net escape 
fractions (i.e. the escape fraction averaged over solid angle) are shown by symbols: yellow (circles), 
blue (squares) and green (triangles) correspond to H I, He i and He ii respectively. Small and 
large symbols corresponding to W^Mq and 10^ Mq Pop HI clusters respectively. For comparison 
black (double dot-dashed) lines show the escape fraction for H I when no helium is included in the 
calculation, with the black circle showing the corresponding angle- averaged escape fraction. 
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Radius; r [kpc] 



Fig. 3. — The escape fraction (thick hnes) and ionization state (thin Hnes) as a function of distance 
from the source ak, = 85.45° (the angle at which the H I escape fraction drops precipitously — see 
Fig. [2|, in the Dove & Shull (1994) Gaussian model for a 10^ M© Pop III cluster. Yellow (sohd), 
blue (dashed) and green (dot-dashed) correspond to H I, He i and He ii respectively. 
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of the cloud population (thin lines) and the average over a much larger number of realizations 
(thick lines). Once again, we use clouds of 30 pc radius, with a volume filling factor of 20% and 
a clumping factor of 10. There is significant variation in escape fraction between realizations — the 
angle-averaged escape fraction for He iii varies from about 55% to 70% with the mean being close 
to 60%. 

To explore the effects of cloud size, the right panel of Figure [5] shows the escape fraction for 
the same model but now using clouds of 3 pc radius. As expected, there is much less variation 
in the escape fraction when using these smaller clouds — since the volume filling factor remains 
the same, the number of clouds is increased by a factor of 1000, reducing the fluctuation from 
realization to realization. The mean escape fractions are not significantly changed relative to the 
case of 30 pc clouds however. Note also that helium is especially sensitive to the presence of clouds, 
with individual realizations departing significantly from the average. 



4.3. Variations in Source Spectrum 



We now consider AGN and Pop III starburst cases in the Dove Sz Shull (1994) Milky Way 
background density profile described earlier, with a cloud filling factor of 20%, cloud radius 30 pc 
and clumping factor of 10. Additionally, we now include cases with and without X-rays (in the 
latter of which the source spectrum is cut off above 120 eV). The AGN spectrum is modeled from 



Hopkins et al. (2007) — this gives the unabsorbed, intrinsic spectral energy distribution (SED) of 
AGN as a function of their bolometric luminosity which we in turn compute from an assumed black 
hole mass, Eddington limit and accretion efficiency (we assume black holes to be radiating at 10% of 
their Eddington luminosity). For this AGN source, we find (Qh i, Qkc i, Que a) = 1.17, 0.68, 0.41 x 
10^° photons s-^ for a 10*^ M© black hole and (Qh i, Que i, Quc a) = 276, 114, 43 x 10^° photons s"^ 



for a 10 black hole. The Pop HI starburst spectrum is taken from Tumlinson et al. (2004), 
which gives (Qh hQue uQuc a) = 38.10,22.97,1.30 x 10^° photons s"^ for a lO^M© burst and 
(Qh I, Que i, Que ii) = 3810, 2297, 130 x 10^° photons s"^ for a W^Mq burst. 

In some calculations, we place these sources in a dark matter halo of mass IO^Mq at z = 10. 
Such a halo will contain of order IO^Mq of baryonic material. If this halo were able to convert 



baryons into stars as efficiently as the most efficient halos at z = (Leauthaud et al. 2012), this 
would imply a total stellar mass of order lO^M©, equal to our more massive starburst model. In 
reality, the efficiency of both black hole and star formation is likely much lower in such a halo (due 
to feedback effects), so we additionally explore W^Mq starburst and AGN models in such a halo. 
Even this is arguably a very high mass for a black hole at such redshifts in a W^Mq halo. As such, 
our derived escape fractions are likely to be optimistic. 

We display plots of the escape fractions as a function of angle for a number of cases that 
highlight key trends when considering a range of source masses, and the presence of X-rays. These 
trends include: 
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Fig. 4. — Left panel: The escape fraction as a function of angle, 9, for a single realization of the Dove 
&: Shull ( 1994) Gaussian model for a 10^ Mq Pop III cluster with clouds of radius 30 pc, filling factor 



20% and clumping factor of 10. Yellow (solid), blue (dashed) and green (dot-dashed) correspond 
to H I, He I and He ii respectively. The net escape fractions (i.e. the escape fraction averaged over 
solid angle) are shown by symbols. Yellow (circles), blue (squares) and green (triangles) correspond 
to H I, He I and He ii respectively. Right panel: The escape fraction (thick lines) and ionization 
state (thin lines) as a function of radial distance from the source at 9 = 85.45° for the same model. 
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Fig. 5. — Left panel: The escape fraction as a function of angle, 9, for ten realizations (thin lines) 
and the average over many realizations (thick lines) of the Dove & Shull (1994) Gaussian model for 
a 10^ M0 Pop HI cluster with clouds of radius 30 pc, filling factor 20% and clumping factor of 10. 
Same notation as Figure |4] Small symbols show mean escape fractions for individual realizations 
while large symbols show the mean escape fraction averaged over a large number of realizations. 
Right panel: The equivalent result for the same model but with clouds of radius 3 pc. 
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The impact of AGN and Pop III sources is quite different, and the effects of switching off 
X-rays in these two cases differ significantly, because the AGN are more X-ray dominated; 



The difference between a blackbody spectrum (as used in ^4.1[ ) and these sources is substantial, 
owing to the Q-values being significantly higher (lower) for the blackbody spectrum relative 
to the AGN (Pop III starburst) - compare, e.g., the Q values for the 10^ Mq AGN or stellar 
cluster with the relevant values in Section 4.1. 



» Once the source mass exceeds about 10^ M0, the background gas (Dove h Shull 1994) "sat- 
urates" and the escape fraction approaches unity for H and He; 

» Interestingly, X-rays boost the escape fractions overall, and are particularly significant for the 
lower-mass AGN or starburst masses. We also see clearly the advancing of the He iii I-front 
relative to that of H 11 or He 11, when X-rays are present. 

In summary, the values of /esc found in the above subsections using the background density 



from Dove & Shull (1994) are quite high, more than 80-90% in some cases for H I and He i. These 
cases are meant to be illustrative rather than exact as we are placing a variety of sources (including 
primordial stars and quasars) in what is essentially a Milky Way-type disk. Some of these sources 



also have significantly higher Q- values than that assumed for a single 0-type star, 10^^ s ^, inlDove 



&: Shull (1994). In the following subsection we will explore primordial sources embedded in dark 



matter halos appropriate to the early Universe. 



4.4. Variations in Galaxy Profile and/or Source Location 



Finally, we consider two alternatives to the Dove &: Shull (1994) galaxy profile: the standard 



NFW profile (Navarro et al. 1997 with the gas profile appropriately scaled to the dark matter 



density through cosmological parameters), and the /3-profile which we use to approximate the 
density distribution of gas in hydrostatic equilibrium at the virial temperature in an NFW halo 



(Ricotti & Shull 2000). For the NFW profile, the density normalization assumes that the halo 



has a virial density contrast as predicted by the spherical top-hat collapse model for our chosen 
cosmology {VLm = 0.275, JIa = 0.725; Komatsu et al. 2011), and the scale radius is set using the 



fitting function of Gao et al. (2008). For the /3-profile, we set parameters using the results of Ricotti 



k, Shull ( 2000 ) such that the profile approximates that of gas in hydrostatic equilibrium in an NFW 



density profile. 

We explored cases for a range of galaxy dark matter masses but display here the results for a 
10^ Mq halo, whose virial radius is about 1.4 (physical) kpc at z = 10. In the Dove Sz Shull (1994) 



model used in earlier sections, the density profile drops away from the source as a Gaussian; thus, 
the escape fraction quickly reaches a constant value with distance from the source. In the case of an 
NFW or /3-profile, however, the density falls off only as r~'^ on large scales, so the escape fraction 
converges more slowly. We therefore use the virial radius of the halo as the radius at which to stop 
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Fig. 6. — Left panel: The escape fraction as a function of angle, 0, averaged over all realizations 
of the Dove & Shull ( 1994 ) Gaussian model for a W^Mq AGN with clouds of radius 30 pc, filling 
factor 20% and clumping factor of 10. Same notation as Figure [4} Right panel: The result for the 
same model but with X-rays removed from the SED. 
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Fig. 7. — Left panel: The escape fraction as a function of angle, 9, averaged over all realizations 
of the Dove & Shull ( 1994 ) Gaussian model for a W^Mq AGN with clouds of radius 30 pc, filling 
factor 20% and clumping factor of 10. Same notation as Figure [4} Right panel: The result for the 
same model but with X-rays removed from the SED. 



-17- 



integrating and measure the escape fraction. The virial radius marks the dividing hne between the 
galactic environment and the IGM. Therefore, any photons escaping beyond the virial radius are 
able to ionize the IGM, in keeping with the usual definition of the escape fraction. 

For these two density profiles (NFW and /3), we consider cases where the source is at the 
center, and cases where the source is placed 100 pc from the halo center (~ 10% of the halo virial 
radius). The figures displayed in this section show the escape fraction at the virial radius as a 
function of angle, as in earlier sections of the paper. For the off-center source cases, note that 
9 = —90° {9 = +90°) is the direction from the source directly through (directly away from) the 
center of the halo, and, consequently, always has the lowest (highest) escape fraction. The angle 
= 0° corresponds to directions perpendicular to the halo center from the source. Most of the 
solid angle is around = 0°, so this typically gives the best estimate of the mean escape fraction. 

Considering first the case of a /3-profile with a lO^M© mass AGN black-hole source (Figure g, 
there is a large difference when X-rays are excluded from the AGN SED, indicating that it is 
primarily the X-rays that are able to penetrate the halo and escape into the IGM. Furthermore, 
in the case including X-rays, it is the He II continuum photons (energies > 4 Ryd) which have the 
highest escape fraction, due to the low cross-section for photons at these energies. For the 10^ M© 



Pop III starburst (Figure 10), /esc is nearly unity for H I and He II ionizing radiation for all of the 
escape directions for a /3-profile. The /esc values for He ii ionizing radiation is significantly lower 
for all angles, approaching zero towards the direction of the galaxy center. If the source is instead 
placed precisely at the halo center, the escape fraction (which is now independent of angle due to 



the spherical symmetry) corresponds to that at = —15° in Figure 10 

In the cases of the 10^ M0 AGN, the values of /esc saturate quickly for all species to 1, except 
for the NFW cases near the very core (owing to the high densities there). X-rays make a difference 
in the direction of center but not otherwise. Hence, we do not display them here, although here 
too, the He curve gets ahead of H in the very core of the galaxy. We find similar results for the 
case of the 10^ Mq Pop HI starburst, except that there is very little difference with X-rays included 
owing to the low X-ray contribution to the SED in this case. 



Figure 11 shows the dependence of escape fraction on radial distance for a IO^Mq AGN source 
placed 100 pc from the center of our /3-profile model. The radial dependence is shown along three 
different directions — directly away from the halo center, perpendicular to the direction to the halo 
center, and directly through the halo center. In the left panel, which includes X-rays, the escape 
fractions quickly asymptote to almost constant values of between 10% and 60% in all directions. 
The right panel of Figure [Tl] shows the same calculation but with X-rays removed from the source 
SED. The behavior is substantially different in this case. In directions through the halo center 
and perpendicular to that direction the escape fractions plummet rapidly, indicating that it was 
the X-rays which were able to escape through these directions most easily. Directly away from the 
center, the escape fractions remain significant, but continue to fall out to the halo's virial radius. 



In summary for this subsection, we find that an NFW gas distribution produces much lower 
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Fig. 8. — Left 'panel: The escape fraction as a function of angle, 0, averaged over all realizations of 
the Dove & ShuU (1994) Gaussian model for a lO^M© Pop III cluster with clouds of radius 30 pc, 
filling factor 20% and clumping factor of 10. Same notation as Figure [4} Right panel: The result 
for the same model but with X-rays removed from the SED. 
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Fig. 9. — Left panel: The escape fraction as a function of angle, 0, averaged over all realizations 
of a /3-profile appropriate to a W^Mq NFW halo at z = 10 for a W^Mq AGN placed 100 pc from 
the halo center and with clouds of radius 30 pc, filling factor 20% and clumping factor of 10. Same 
notation as Figure |4j Right panel: The result for the same model but with X-rays removed from 
the SED. 
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Fig. 10. — Left panel: The escape fraction as a function of angle, 9, averaged over all realizations of 
a /3-profile appropriate to a IO^Mq NFW halo at z = 10 for a lO^M© Pop III cluster placed 100 pc 
from the halo center and with clouds of radius 30 pc, filling factor 20% and clumping factor of 10. 
Same notation as Figure |4j Right panel: The result for the same model but using the NFW itself 
rather than the corresponding /3-profile. 
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Fig. 11. — Left panel: The escape fraction (thick lines) and ionization fraction (thin lines) as a 
function of distance, r, from the source, averaged over all realizations of a /3-profile appropriate to 
a IO^Mq NFW halo at z = 10 for a IO^Mq AGN placed 100 pc from the halo center and with 
clouds of radius 30 pc, filling factor 20% and clumping factor of 10. Same notation as Figure |4j 
Right panel: The result for the same model but with X-rays removed from the SED. 
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escape fractions overall for both H and He relative to the /3-profile, particularly through the halo 
center where the density rises rapidly. Additionally, the presence of X-rays boosts /esc for H I, He i 
and He il; this is especially dramatic for He ii. 



5. CONCLUSIONS AND DISCUSSION 



We have calculated the escape fraction of ionizing radiation, /esc, for both hydrogen and helium 
using ID semi-analytic models of primordial galaxies. We have considered many scenarios that vary 
the galaxy density profile, the source type (first stars or AGN), location, and spectral index, as 
well as gas clumping factors. Key new points in our work include calculations of /esc for He i 
and He ii ionizing radiation, and the inclusion of the process of partial gas ionization by X-rays. 
We find that X-rays play an important role in boosting the /esc values for H I, He i and He ii, 
especially He ii. Escape fractions from primordial halos containing a bright AGN may reach as 
high as 30:20:10% for H ii:He ii:He iii (mostly due to the significant contribution from X-rays in 
such sources), while a burst of Pop. HI star formation could result in escape fractions as high as 
90:90:40% for H ll:He ll:He III. For sufficiently hard first-light sources, the helium ionization fronts 
closely track or even advance beyond that of hydrogen. 

These calculations have an impact on a number of different observational probes of first-light 
sources and reionization models. First, the current WMAP7 results indicate that the Thomson 
optical depth, Te ~ 0.088, allowing for an extended period of partial ionization of H and /or He 
beyond z ~ 6. As discussed in Shull & Venkatesan (2008), the need for an extreme reionization 
scenario is reduced for H, if /esc,He is competitive with /esc,H- With primordial helium at 1" = 0.25 
by mass, the number ratio y = nue/'^H = 1/12. The inclusion of helium in models typically 



increases Te by 16% (Shull &: Venkatesan 2008), accounting for singly ionized helium at 3 < z < 7 



and doubly ionized helium at z < 3. 

Second, the values of /esc and its redshift evolution can impact the evolution of the optical 
depth of He ii and the ratio of He iii to He ii at z ~ 2.5-3 (iKhaire & Srianandl 120121 IShull 



et al. 2012b), during late He ii reionization. The ratio of He II to H I, denoted by t], has been 
observed to have factor-of-ten variations on Mpc scales, along the sightlines to several "He-II 
quasars" at z ~ 2.4-2.9 ( |Shull et al.||2004[ |2010t [Fechner et ar]|2006| [Syphers et al.||20"T2l ). Because 



He II is photoionized by 4-Ryd continuum radiation, these ?7-variations are probably dominated by 
fluctuations in the intrinsic spectra of the AGN sources, modulated by He II opacity changes in the 
IGM (primarily from Lyman-limit systems) . Ultraviolet spectroscopic studies of low-redshift AGN 
(Zheng et al. 1997; Shull et al. 2012b) find rest-frame (1-2 Ryd) ionizing spectra, F^, oc with 
a wide range of spectral indices from hard (a ~ 0.5) to soft (a ~ 3). This is consistent with the 
critical spectral index derived in Section 3 as well as in the Appendix below when accounting for 
frequency-dependent effects in the source spectrum. The composite EUV spectra of low-redshift 
AGN have recently been measured to have (a) ~ 1.4 — 1.5 between 1.0-1.5 Ryd ( Shull et al.|2012"b ). 
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Last, calculations of /esc impact the converse problem of radiation trapped in primordial galax- 
ies, and have consequences for detecting these galaxies through emission- line signatures. Such 



emission-line signatures of very hot, low metallicity massive stars (Tumlinson &; Shull 2000 Tum- 



linson et al.||200T ) include strong He II recombination lines such as A1640 and A4686. 



Looking ahead, science drivers for future UV spectroscopic missions (Tumlinson et al. 2012 



McCandliss et al.|2012 ) include telescopes and spectrographs with spectral coverage down to the hy- 
drogen Lyman limit (912 A). This coverage will allow one to directly measure escaping (hydrogen) 
LyC radiation from low-redshift starburst galaxies. This far-UV band also probes He II Gunn- 
Peterson absorption between z ~ 2 — 3, an interval that spans the patchy He II post-reionization 



epoch in the rest-frame of the intervening IGM (Shull et al. 2004 2010 Syphers & Shull 2013) 



Observations of "He li quasars" at somewhat higher redshifts, z = 3.0 — 3.5 ( Syphers et al.|[2009b|a 



2012 Worseck et al. 2011) can observe the He II absorption during the epochs when the He III 



ionized bubbles have not yet overlapped. Placing direct limits on /esc,He faces significant observa- 
tional challenges, however, and the likelihood of a FUV instrument as a successor to COS remains 
uncertain at present. 
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populations, and Yue Zhao, Yang Wang, and Long Yan Yung for useful input. AV gratefully 
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A. Analytic Calculations of Escape Fractions 

In order to derive the values of /esc for a general spherically-symmetric density profile, we 



begin with Equation (12) of|Dove fc SImlll (|1994D, 



/ei 



1 



47raB 
5*0 



nennr^dr . 



(Al) 



This preprint was prepared witii tlie AAS lAT^rjX macros v5.2. 



Fig. 12. — The coordinate system within a galaxy halo assumed for the analytic derivations in the 
Appendix. The light source (AGN or star cluster) is at a distance D from the galaxy center, R is 
the distance from the center to the current position, and r is the position difference between D and 
R. Note that 9 + (j) = 90°, and that 9c is the critical opening angle of the cone of escaping ionizing 
radiation, as defined in earlier sections. 
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Although we eventually want to replace Sq with a frequency-dependent function for the source 
luminosity, for now we will assume it is a constant. For a general spherically-symmetric galaxy 
density profile, the integral in the second term above reduces to ( Ricotti Sz Shull|[2000 ) : 



nc,c"-H,c^^dr 

l + [D'^+r^ + 2Dr s\n{(l))]/Rl ' 



(A2) 



where, as shown in Figure 12 D is the distance of the source from the galaxy center, R is the 
distance from the center to the current position, and r is the position difference between D and R 
(with R^ = D'^ + r'^ + 2Drsm[(p]). Note that, as shown in Figure [T2I 6 = 90° -(j). The subscript "c" 
in this section represents the core values of each quantity. We assume that the baryonic distribution 
is given by the /3-profile frequently used to model the gas distribution in galaxies (and used in earlier 
sections of this paper), and does not lead to a central cusp in the density (as in the case of an NFW 
profile) : 



and the electron density is given by, nc{R) 
doubly ionized. 



3P 

2 



l + WRcr\ ' ^^^^ 

nii{R) + 2nHe(-R) oc nii{R), assuming that helium is 



We consider the source (stars or QSO) to be at the center of the halo, so that D 
r = R, and also take 3/3/2 ~ 1. Then the integral above reduces to one purely in r. 



and 



where a = 1, and b = R^. 



r'^dr 



'0 l + ii-J' 
This has the solution: 
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(A4) 



(A5) 



If we evaluate this at the center (r = 0) and at the galaxy halo's edge (r = R^), the general 
analytic expressions for the escape fractions in a /3-profile for H and He without source frequency 
dependence, and with the source at the halo center, are respectively: 



esc,H 



and 



esc, He 
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47raB,H 

'S'o.H 

47raB,He 
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(A6) 



(A7) 



Note that He in this section always refers to the double ionization of helium. For a /3-profile, the 
ratio R\i/Rc is essentially a concentration parameter. The second term on the right-hand side in 
each of the two above expressions has several common factors, but this second term will in general 
be smaller for the case with H (see discussion of VFh vs. Wne in Section 4.1), implying that for 
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most sources, /esc will be greater for H than for He. Note that this does not account for the critical 
role of X-rays in the early stages of helium and hydrogen ionization. 

We compare the values implied by the above expressions with our numerical results for a 
central source in a /3-profile; in contrast to the numerical results shown in Section 4, we exclude 
X-rays and clumping in this comparison so as to match the conditions in Equations (A6) and (A7). 
In these cases, Rc = 0.28 kpc and the concentration parameter is about 4.9. We find that the 
analytic expressions reproduce the numerical results to within a few percent for the higher-mass 
sources like the 10^ M© AGN or Pop III star cluster. However, the analytic expressions for /esc 
match some of the values for the lower-mass (10^ M©) AGN or star cluster poorly, sometimes even 
leading to negative values, particularly when the numerical results yield very low /esc values close 
to zero. This discrepancy likely arises from not including X-ray ionization in the analytic formulae, 
an effect which becomes important at ionization levels below ~ 10%. The analytic expressions are 
also very sensitive to the exact balance of number densities between the various H and He species, 
relative to the ionizing photon fluxes. 

One possible way to include clouds in the analytic formalism is by noting that the addition of 
clouds does not change the total column density of atoms and ions along the line of sight, but that 
the recombination timescale, going as ~ l/(ane), will differ between cloud and non-cloud regions. 
The second term on the RHS of Equation (Al), Wj for each species as defined earlier in Section 4.1, 
goes as the product ngrij; within cloud regions, Ug and ni are increased by a factor (l/[/ + (1 - 
f)/C]), while between clouds He and rii are decreased by a factor (l/[/(C — 1) -|- 1]). Recall from 
earlier sections that / is the cloud volume filling factor and C is the cloud overdensity relative 
to the non-cloud gas. Thus, clumping factors could be reproduced in the analytic expressions by 
multiplying Wi by the factor: //[/ + (1 - f)/C]^ + (1 - /)/[/(C - 1) + 1]^. 

For our typical choices in the paper (/ = 0.2, C = 10), this factor is approximately 2.65. 
This improves some of the discrepancies between the analytic and numerical results for the low- 
luminosity sources mentioned above, but not to the desirable range of being within a few to 10 
percent of each other. 

Now, wc assess the impact of frequency dependence in the input spectrum, and replace Sq 
with So{u/uo)~^"~^^\ Thus, the first part of the second term in Equation (Al) becomes: 

47raB 
JduSo{u/uo)-(<^+^) 

The constant So has to be evaluated through integration from a threshold energy, e.g., from huH = 
1 Ryd or hi/ue = 4 Ryd, up to infinite energy, and by setting the integral equal to the H-ionizing 
(or He-ionizing) photon flux in units of photons s~^. If we evaluate the above integral between the 
ionization threshold energy of interest up to energies hvi = A^huo), where vi is a constant multiple 
of vo, then the integral reduces to (switching variables to vjvo): 
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uoSoil - A-'^], forH. 



(A9) 



uoSo (^) (^) ^ ^ ^ = z^o5o[4-° - A-% for He . (AlO) 

For H I, huQ = 1 Ryd, whereas for He ii, hvQ = 4 Ryd. Note that the constant 5*0 through its 
normahzatfon is always proportional to the product aSo^u/Woi^ — A~°')] for H, and to the product 
a5o,He/ [2^0(4"" — for He. When integrating up to infinity (i.e., A = 00), the last term on 

the right-hand side simply goes to zero. Therefore, the general analytic expressions for the escape 
fractions in a /S-profile for H and He with source frequency dependence, and with the source at the 
halo center, are respectively: 
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To compare with the calculation of the critical spectral index of a 



Rc 



(All) 
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1.8 derived earlier in Section 
3 when H and He photoionizations became equal, we set the second terms on the RHS of the two 
above equations equal to each other. This leads to: 



0.08 I ) , or Oeff - 0.5, 



(A13) 



for the standard values for the recombination coefficients. Here, a^s is the critical source spectral 
index after gas filtering and transmission effects are taken into account. This value of OeS ~ 0.5 is 
lower than the Ocrit ~ 1-82 derived in Section 3 which did not account for these effects. Both these 



values, however, encompass the range found in UV spectroscopic studies of low-z AGN (Zheng 



et al. 


1997; 


Shull et al. 


2012b; 



